Phase Transition in a Random Minima Model: 
g ; Mean Field Theory and Exact Solution on the 
^ ■ Bethe Lattice 



Peter Sollich^, Satya N. Majumdar^, Alan J. Bray^ 

OO ■ 1 Department of Mathematics, King's College London, London WC2R 2LS, UK 

^ Laboratoire de Physique Theorique et Modeles Statistiques, Universite de 
Paris-Sud, 91405 Orsay Cedex, France 

^ School of Physics and Astronomy, University of Manchester, Manchester M13 9PL, 
UK 






C/3 



cd 



O 
o 



o 



E-mail: peter . sollichOkcl .ac.uk 



Abstract. We consider the number and distribution of minima in random landscapes 
defined on non-Euclidean lattices. Using an ensemble where random landscapes are 
'"^ . reweighted by a fugacity factor z for each minimum they contain, we construct first 

G I a 'two-box' mean field theory. This exhibits an ordering phase transition at Zc = 2 

above which one box contains an extensive number of minima. The onset of order 
is governed by an unusual order parameter exponent /3 = 1, motivating us to study 
the same model on the Bethe lattice. Here we find from an exact solution that for 
^ ' any connectivity fi + 1 > 2 there is an ordering transition with a conventional mean 

^^D , field order parameter exponent (3 — 1/2, but with the region where this behaviour is 

rJ: ' observable shrinking in size as l//i in the mean field limit of large /i. We show that 

^b . the behaviour in the transition region can also be understood directly within a mean 

1^ I field approach, by making the assignment of minima 'soft'. Finally we demonstrate, 

f— ^ ' in the simplest mean field case, how the analysis can be generalized to include both 

OO , maxima and minima. In this case an additional first order phase transition appears, 

to a landscape in which essentially all sites are either minima or maxima. 



. F^. , 1. Introduction 



The statistics of the number of stationary points (maxima, minima and saddles) 
in a random landscape plays an important role in understanding both the static 
and the dynamical properties of many systems such as structural glasses [1], spin 
glasses [2], clusters and biomolecules [3], continuum percolation [4], rugged landscapes 
in evolutionary biology [5], quantum cosmology [6], string theory [7, 8], particles 
in a random potential [9, 10], and also in several associated problems in random 
matrices [8, 11, 12, 13]. In particular, the statistics of the total number of minima 
is important to understand in the context of glassy materials where the system typically 
gets trapped for a long time in a local minimum of the energy (or free energy) 
landscape [1, 2]. 
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In theoretical studies one typically models a random energy landcsape as a 
smooth random manifold (typically Gaussian) sitting on an underlying continuous 
space. For such smooth Gaussian random surfaces in a continuum, there is a finite 
density of local minima (expected number of local minima per unit volume) which 
can be computed from the celebrated Kac-Rice formula [14] and its multi-dimensional 
generalizations [15, 16, 17]. Recently these formulae have been extended to compute 
the expected number of saddle points at a given fixed energy and also with a fixed 
index number of saddles [18, 19, 20]. Similarly the variance and higher moments can 
also be computed in principle. However, this machinery is not easily extendable to 
cases (i) where the surface is non-Gaussian and (ii) where the underlying space on 
which the energy landscape resides is discrete, e.g. a regular Euclidean lattice. The 
latter case is particularly relevant in the practical context of numerical simulations of 
the random energy landscape where one is obliged to discretize the underlying space. 
Hence it is interesting to study the statistics of the number of local minima in lattice 
models of energy landscapes, in particular where the energy distribution at each site is 
non- Gaussian in general. 

With these two motivations in mind, a simple lattice model of an energy landscape 
has recently been introduced (hereafter referred to as the 'random minima' model) [21]. 
In this random minima model, a random energy Ei sits at site i of a lattice of A^ sites 
with periodic boundary conditions. The energies Ei are drawn, independently from site 
to site, from a common continuous distribution p{E), not necessarily Gaussian. Any 
such choice of the set {Ei} defines a realization of a random energy landscape. For a 
given realization, a site i is a local minimum if Ei < Ej for all sites j which are nearest 
neighbours of site i. Let M denote the total number of local minima. Evidently M 
will vary from one realization of the landscape to another and one is interested in the 
probability distribution P{M, N) of M for a given size A^ of the system. The energies 
at different sites are uncorrelated in the random minima model. Hence it can be viewed 
as an effective 'coarse grained' lattice model of a continuous random energy manifold 
in the limiting situation where the correlation length between the energies at different 
points in space is smaller than the lattice spacing. 

The advantage of this simple model is that many questions regarding the statistics 
of the number of minima are analytically tractable [21]. Besides, the same distribution 
P{M, N) has appeared recently in seemingly unrelated problems such as random 
permutations [22, 23], ballistic deposition models [24] and also simple models of 
glasses [25]. The distribution P{M,N) turns out to be strictly universal, in the sense 
of being (even for finite A^) independent of the on-site energy distribution p{E) as long 
as the latter is continuous [21, 24]. This is not difficult to see: transform from the Ei 
to new variables Xi = q{Ei), with q{E) = J_^dE'p{E') the cumulative distribution 
function. For continuous p{E) , this transformation is monotonic so that minima in the 
ii^j- landscape are identical with minima in the landscape defined by the Xj. But the Xi 
are uniformly distributed in the interval [0, 1], so it suffices to consider this particular 
distribution - denoted Q{x) below - to obtain P{M, N) for any continuous p{E). 
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The average number of minima can be trivially computed, (M) = N/[{fi + 1) + 1] 
where /x + 1 is the co-ordination number of the lattice. For example, for a rf-dimensional 
hypercubic lattice, n = 2d — 1 whereas for a Bethe lattice fi is just the branching ratio. 
Similarly, the variance of M can also be computed exactly for various lattices such as a 
1-d chain [21, 24], the 2-d square lattice [24] and the Bethe lattice [21]. The distribution 
P{M, N) has a Gaussian peak near its mean (of width ~ vN), but a non-Gaussian tail 
far from the mean. The non-Gaussian tail is described by a large deviation function 
that can be computed exactly in 1-d [21]. Also, on any given lattice, M can at most 
take a value Mmax- This follows from the fact that if a site is a local minimum, none of 
its neighbours can be a local minimum (nearest neighbour minima exclusion principle). 
For example, on a bipartite lattice, consisting of two 'boxes' each containing A^ sites 
and where each site has nearest neighbour connection to all sites in the other box, one 
cannot have minima in both boxes and so Mmax = N. In Ref. [21], the probability of 
the maximal packing configuration P(Mmax5 -A) was studied and was shown to decay 
for large A^ as P{M^a_^,N) ~ 7~^, where the constant 7 was exactly computed for a 
number of lattices. 

The purpose of this paper is to go beyond the 'counting problem' in the random 
minima model and study its thermodynamics and the associated phase transition by 
introducing a fugacity z for each local minimum. For this purpose, the relevant object 
of interest is the generating function (or grand partition function) 

G{z,N) = Y,P{M,N)z'^ (1) 

M 

and the associated equation of state. The latter tells us how the density of minima 
p = {M)/N depends on z, where the average is over the original ensemble of random 
landscapes reweighted by a factor z^ for each configuration. 

At this point it is important to note that due to the nearest neighbour minima 
exclusion principle, the thermodynamics of the random minima model is similar 
in spirit, though not in its details, to the well studied 'hard sphere lattice gas' 
model [26, 27, 28, 29, 30, 31]. In the latter model, when a molecule or a hard sphere 
occupies a lattice site, a similar exclusion principle holds in that all neighbouring 
sites have to be empty. If W{M, N) denotes the number of ways of putting M hard 
particles (with this constraint of nearest neighbour exclusion) on a lattice of A^ sites, 
the corresponding grand partition function is defined as 

Z{z,N) = J2^{M,N)z^. (2) 

M 

Note the important difference between the two models. In the hard sphere model with 
M particles, one attaches a uniform weight 1 to each allowed configuration of the M 
particles. On the other hand, in the random minima model, for each configuration 
of M local minima, the associated weight comes from an entropic factor obtained by 
integrating over all possible -Ej's associated with the given configuration of the M local 
minima. 
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The hard sphere model is well known to undergo a thermodynamic phase 
transition as one increases the fugacity z through a critical value Zc in two or higher 
dimensions [28, 31]. For z < Zc, the system is in a low density 'disordered' or 'fluid' 
phase and for z > z^ it is in a high density 'ordered' or 'crystalline' phase. Based on 
the qualitative analogy between the two models one therefore expects a similar phase 
transition from a disordered to an ordered state in the random minima model, also 
upon increasing the fugacity z. Indeed, recent numerical studies by Derrida for a 2- 
d random minima model indicate the presence of such a phase transition [32]. It is 
important to understand whether this phase transition in the random minima model is 
similar/different from that of the hard sphere lattice gas model. 

With this in mind, we study the thermodynamics of the random minima model 
on the Bethe lattice. The hard sphere model was solved exactly on the Bethe lattice 
many years back [31] and has been revisited recently [33]. In this paper we present an 
exact solution of the random minima model on the Bethe lattice which turns out to 
be technically somewhat harder than the hard sphere solution on the same lattice. In 
addition, we study analytically a rather simple mean field theory of the random minima 
model which also exhibits a phase transition at a critical value Zc = 2. We show that in 
the low density phase (for z < 2) the average number of minima is of order unity in the 
thermodynamic limit, while in the high density phase (for z > 2) the average number 
of minima is extensive with a finite density. 

We survey the hard sphere lattice gas briefly in Sec. 2. There is no non-trivial mean 
field theory for hard particles, so we start directly with the Bethe lattice case (Sec. 2.1). 
A mean field theory can be constructed if particles are made soft, i.e. if occupation of 
neighbouring sites is permitted subject to some penalty. As we show in Sec. 2.2, when 
the penalty parameter is made large this approach nicely retrieves the results for the 
Bethe lattice in the limit of large connectivity. 

We turn to our main subject, the random minima problem, in Sec. 3. Here there 
is a non-trivial two-box mean field theory and we discuss this first, in Sec. 3.1, and also 
extend it to study the joint statistics of the number of minima and maxima (Sec. 3.2). 
Next we analyse the random minima problem on the Bethe lattice (Sec. 3.3) and finally 
we consider a mean field theory with soft assignments of minima in Sec. 3.4. Again 
we will see that these two approaches give the same results in the respective limits 
of large connectivity and almost-hard assignments. Considering these limit cases also 
helps to clarify why the direct mean field approach of Sec. 3.1 gives an unusual apparent 
order parameter exponent near the phase transition. We summarize and list some open 
questions in Sec. 4. 

2. Hard sphere lattice gas 

In this section we revisit briefly the hard sphere lattice gas model, on a Bethe 
lattice [31, 33] and in a two-box mean field theory with soft particles. This will serve 
to introduce the techniques we will deploy later for the problem of minima in random 
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landscapes. Note that the simplest mean field theory, a fully connected lattice, makes 
no sense as the presence of a single particle would exclude particles from all other sites. 
Also the simplest improvement over this, a fully connected bipartite lattice, is trivial: 
one of the two boxes, i.e. partitions of the graphs, is always empty, and in the other the 
particles are then non-interacting. In the random minima problem, on the other hand, 
already this approach produces a phase transition as discussed in Sec. 3.1 below. 

2.1. Bethe lattice 

Consider first a Cayley tree| with branching ratio /i, of depth /, i.e. with / layers below 
the single root node. Call Zq { the grand partition function constrained to run over 
all configurations that do not (or do, respectively) have a particle at the root. The full 
partition function is them Z^'''^ = Zq +Z} . We do not write explicitly the dependence on 
z, while the superscript (/) indicates indirectly the number of sites N = (yu'+^ — l)/(yU — 1) 
in the tree. 

The quantities Zq { obey the following recursions over the tree depth: 

4'+i) = (4') + zff (3) 

Zf+'^ = z{zPY (4) 

with Zq = I, Zi = z. For example, if no particle is present at the root node of a 
tree of depth / + 1, then the fi sites in the next level of the tree are each allowed to 
be either occupied or not; the partition sum is then the product of fi unconstrained 
partition functions Zq + Z^ for each of the subtrees of depth /. This gives the first 
equation above. For the second equation, one notes that if a particle is present at the 
root then each of the fi sites below must be empty. The partition sum is then a product 
of the appropriate constrained partition sums Zq for the subtrees, with an extra factor 
z to account for the particle at the root. 

The two recursions can be combined into one for the ratio S^''^ = Z\ /Zq , giving 

^^'^'^ = (1 + si^Y ^^^ 

with S**-"-* = z. From S''-^-' one can determine the density at the centre of a Bethe lattice 
of depth / + 1, obtained by connecting /i + 1 Cayley trees of depth / to a central node 
(see Fig. 1). Taking the appropriate ratio of the partition sum with the central site 
occupied to the total partition sum yields 






(7) 



(l + 5W)M+i + ^ ■ 

I We note that in some of the literature what we caU a Cayley tree is termed "rooted Cayley tree" , 
while the term "Cayley tree" is then used for what we call a Bethe lattice, i.e. a tree where every node 
except those on the boundary has /i + 1 neighbours. 
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Figure 1. Sketch of a Bethe lattice witli jj, — 3: every interior node has n + 1 — 4 
neighbours. This lattice can be obtained by connecting /i + 1 = 4 Cayley trees - 
indicated by the dashed lines - of branching ratio /i = 3 (and in this case depth ^ = 1) 
to the central node. 

For low z the recursion for S**^'-* has a single fixed point. The bifurcation to an ordered 
state occurs when dS^'"'^^^ /dS'^^^ = —1 at the fixed point§. This requirement together 
with the fixed point condition itself gives 

Sc = ^, z, = fi^/ifi-ir+\ pc = ^. (8) 

/i — 1 /i + 1 

The divergence oi Zc at fi = 1 makes sense: for /i = 1 we have a chain, which as a one- 
dimensional system with only short-range interactions cannot exhibit a phase transition. 
For z > Zc the recursion for S''-'-' converges to a period-two sequence, S^^^^ -^ S, 
5{2fc+i) ^ s, where 

^ ^ {1+sY' ^ ^ (1 + sr ■ ^^^ 

As anticipated, this means the system is ordered, with alternating layers of the lattice 
preferentially occupied/empty; the densities are: 

z 1 

^ " (1 + Sy+^ + z^ ^i/^5-(^+i)/^ + 1 ^^^> 

_ z 1 

^ " (1 + SY+^ + z^ ^i/i^s-ip+i)/^ + 1 • ^1^) 

(Given the initial condition 5**^°^ = z > Sc^ the even layers should be the occupied ones, 
i.e. S > S.) Mathematically, further bifurcations could occur for larger z, but physically 
this is implausible. In general, the equations for S and S need to be solved numerically. 
For large /i simplifications occur, however. It is to this mean field limit that we now 
turn. 

§ If we write the recursion as 5''''+^^ = f{S^^^), then the bifurcation is to a cycle of two solutions, 
S — f{S) and S — f{S). Near the bifurcation S and S are close, so one can expand 5* = 
f{S) + iS- S)f'{S) + 0{{S - 5)2) which gives {S -S)[l + f'{S)] = 0{{S - S)^) and hence f{S) = -1 
at the bifurcation itself. 



Phase Transition in a Random Minima Model 7 

2.1.1. Large fi, above the transition Here we take fixed z > z^. It is then not hard to 
see that for yU ^ oo one gets S = z, S = z{l + z)~^. (This is self-consistent, since S is 
exponentially small and so (1 + S)^ ^1-) The resulting densities are, to leading order, 

p = l/{z-' + l) = z/{z + l), p = 2(l + 2)-(^+i) . (12) 

This is plausible: p is very small so that the odd layers are basically empty. Then p is 
just determined by the activity z, which attributes weights z and 1, respectively, to the 
configuration with or without a particle on a site of the even sublattice. 

2.1.2. Large p, around the transition For p — > cxd, the critical activity from (8) becomes 
Zc = e/p. We therefore set z = z/p to explore the region around the ordering transition. 
Since the critical density pc = ^/{p + 1) and partition sum ratio ^c = l/(/x — 1) are also 
C(l//i), we put likewise p = p/p and S = S/p. The fixed point in the disordered phase 
then obeys, from the laige-p limit of the fixed point of (5), 

S = ze-^ (13) 

From (7) the density becomes p = ze^^ = S, so the 'equation of state' is simply 

z = pe^ . (14) 

In the ordered phase, on the other hand, one has from (9) 

S = ze-^, S = ze-^ (15) 

with again p = ze~^ = S and p = S. So the activity and the densities in the even/odd 
layers are related by 

2 = peP= ~peP . (16) 

The two densities obey 

pe'~P = ~pe-'p (17) 

and hence the critical point is at p = p = pc = 1, 5c = e as expected. For higher 5, 
the densities deviate from each other with a standard square root singularity, to leading 
order, p — 1 = 1 — p~ {z — z^)^!'^ . This corresponds to an order parameter critical 
exponent /9 = 1/2 as expected for a mean field model. For general z > z^, the last 
two equations - which together determine the equation of state of the ordered phase - 
need to be solved numerically. (One could choose, say, p > 1, find the corresponding p 
from (17), then determine z.) The asymptotic behaviour for z^lisp = z, p = ze~^ 
which matches with the z <^ 1 limit of (12) as it should. 

2.2. Mean field theory with soft particles 

The large connecitivity limit p ^ oo discussed above must correspond to a mean field 
theory that one ought to be able to construct directly, without having to first solve for 
lattices of finite connectivity. As explained above, a fully connected lattice makes no 
sense as the presence of a single particle would exclude particles from all other sites. One 
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is therefore led to considering a fully connected bipartite lattice. This can be thought of 
as two boxes ('left' and 'right') with N sites each; every site is connected to all others in 
the other box. If we now directly enforce the hard repulsion of particles on neighbouring 
(connected) sites, the model is trivial: as soon as one box contains any particles, the 
other one must be completely empty. The density in the non-empty box is then just 
p = z/{z + 1) as determined by the fugacity, and the system is always ordered. 

To retrieve the ordering phase transition, one needs to introduce a soft repulsion. 
Here we give a configuration with M and M particles in the two boxes weight 
z^+^ exp{—aMM/N). Sending a ^ oo then recovers the hard repulsion. 

The partition function for this soft repulsion model is 

Zi^^N) =T.(V)(l) z-^^^^e'-^^'^ (18) 

M,M \ ' \ ' 

and can be evaluated by introducing the densities p = M/N, p = M/N and evaluating 
using steepest descents for A^ -^ oo: 

A^-^ InZ = ma_x {H (p) +H{p) + {p + p)lnz - app] (19) 
with the entropy 

n{p) = -p\np-{l-p)Hl-p). (20) 
The resulting saddle point conditions are 

ln[(l-p)/p]+lnz-ap = (21) 

ln[(l-p)/p]+ln2-ap = . (22) 

We will be interested in the large a limit where the repulsion is 'nearly hard'; when the 
symmetry between boxes is broken, we assume without loss of generality that it is the 
left box that has the higher density, i.e. p > p. 

2.2.1. Large a, above the transition Taking a large at fixed z, we see that to satisfy 
the second saddle point equation to 0{a) one needs p = 2;exp(— ap) to leading order. 
The a-dependent term in the first saddle point equation then becomes negligible, so 
that (1 — p)/ p = 1/z or p = z/{z + 1). This is just the result (12) on the Bethe lattice 
for large p, as expected, and is consistent with the simple expression obtained from the 
balance of the weights of unoccupied and occupied configurations (see above). Note that 
the density of the almost empty box is p = zexp[—az/{z + 1)] to leading order; this 
does not match with the Bethe lattice result if one naively identifies a with p. So only 
the leading order densities {z/{z + 1) and 0) match while the subleading (exponentially 
small, in the nearly empty box) corrections are not related. 

2.2.2. Large a, around the transition Here we set p = p/a and p = p/a, by analogy 
with the large p treatment on the Bethe lattice. This gives for a ^ oo the saddle point 
equations 

ln{z/p) -~p = 0, Iniz/p) - p = . (23) 
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These can be rewritten as 

z = pe'p= ~peP (24) 

which is exactly the same as the large-/i result on the Bethe lattice. So near their 
respective transitions the two models behave identically, demonstrating that the fully 
connected two-box model with soft repulsion captures the same physics as the Bethe 
lattice for high connectivity. 

3. Random minima 

In this section we turn to our main subject, the arrangement of the local minima of a 
random function on a Bethe lattice. As explained in the introduction, we can without 
loss of generality take the function value at each site i to be a random variable Xi sampled 
from a uniform distribution Q{x) over [0, 1]. We will define binary indicator variables 
nii, setting rrij = 1 if site i is a minimum, i.e. if none of its neighbours has a larger x; 
otherwise we set rrij = 0. (With this convention, leaves of a tree are counted as minima 
if their x is smaller than that of the parent node directly above.) The mi are analogous 
to hard particle occupation numbers since no two neighbouring nodes can have m = 1; 
but the random values Xi introduce other, non-trivial correlations. As in the hard sphere 
model we will multiply the weight of any configuration of the mi, produced by a random 
draw of the Xj, by a fugacity factor z'^ , where now M = J^i^^i is the total number of 
minima. We wish to calculate the density of minima as a function of z, and understand 
whether an ordering transition does again take place for sufficiently large z. 

3.1. Mean field theory 

We begin with the simplest calculation, which is the two-box mean field theory. Each 
box contains A^ sites as before, and the particles in each box are regarded as neighbors 
of all the particles in the other box. Clearly, all sites that are minima must belong to 
the same box. Call the random variables in the left box Xi and those in the right box 
Xi, with i = 1, . . . ,N . 

To work out the generating function (1) we need to find P{M, N), the probability 
of having M minima in our system. Assume first that the minima are in the left box. 
There are M minima if precisely M among the Xi are smaller than all of the Xj, i.e. 
smaller than x_, the smallest of the Xj. Given that the Xj are uniformly distributed 
over [0, 1], the probability distribution of X- is 

P{x^) = N{l-x^f-\ (25) 

so the probability that M of the ojj's are smaller than a;_ is 

where the average is taken using the distribution P{x-). Accounting for the 
configurations where the roles of the two boxes are swapped, the probability of having 
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M minima is twice as large: 

P{M, N) = 2N (^) I' dx_ x5(l - x.f''-^-' (27) 

^ iV!(2iV-M-l)! 

{N - M)\{2N - 1)\ ■ ^ ' 

Note that the value M = is impossible as there is always at least one minimum present, 

and accordingly one has the normalization Y,m=i P{^i ^) = 1 as is easily checked. 

There is in fact a simple counting argument that leads directly to the result (28). 

To generate the A^ random numbers in each box, we can first sample 2N random 

numbers yi, ■ ■ ■ ,y2N, again from the uniform distribution over [0, 1]. We then take a 

bag containing A^ labels 'left' and A^ labels 'right' and, for each of the y^, pull out one 

label from the bag and put yi in the relevant box. Because the order in which we consider 

the different y^ for labelling is irrelevant, we can in particular take them to be ordered, 

1/1 < ... < y2N- Then a configuration with M minima in the left box is one where yi to 

yM have got labels 'left' and yM+i the label 'right'. Keeping track how many 'left' and 

'right' labels remain in the bag at each step of the labelling, and including the overall 

factor of 2 for the reverse situation where the M minima are in the right box gives 

^/,. ,.x N N-1 N-(M-l) N 

P(M,N) = 2 \ '- ■ 29 

^ ' ^ 2N 2N-1 2N-{M-1) 2N - M ^ ' 

m {2N-M)\ N 

{N-M)\ {2N)\ 2N-M ^ ' 

(2N -M -l\ (2NY^ , , 

= H A'-l )(iv) ■ Pl' 

The second expression is the one most easily seen to agree with (28). The third one 
gives another way of thinking about the result: having M minima (in the left box) fixes 
the first M + 1 labels, and the probability is then the number of arrangements of the 
remaining labels divided by the number of arrangements of all labels. 

Returning now to our original aim of computing the generating function G(z,N), 
it is in fact most convenient to employ the integral form (27) of P{M,N). Then the 
sum over M in Eq. (1) can be evaluated to give 

G{z, N) = 2N r dx_ (1 - x_f-^ {[1 + (^ - l)x-f - (1 - x_f] . (32) 

The following analysis shows that there is a phase transition, at Zc = 2, between a phase 
where the number of minima remains finite as A^ ^ oo and a phase where the minima 
are extensive in number. We begin by setting X- = v/N in Eq. (32), and taking the 
limit A^ ^ oo at fixed v. In this limit the integrand becomes exp[— f (2 — z)] — exp(— 2f ) 
and the upper limit on v tends to infinity, to give 

G(^,Ar^oo) = ^-i = ^, z<2. (33) 

The same expression can be obtained by noting that, for finite M and large A^, the 
number of minima has the geometric distribution P{M, N —^ oo) = 2~^. The result (33) 
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diverges aX z = 2. This indicates that the underlying assumption - that values of x_ 
of order 1/A^ (corresponding to M of order unity) dominate the integral in (32) - no 
longer holds, and suggests a phase transition at Zc = 2. 

For z > 2, the integral can be evaluated using the method of steepest descents. For 
this purpose we write the integral (32) in the form 

G{z,N) = N [^ ^^—{{l-x_)[l + {z-l)x_]}'' . (34) 



/o 1 — a;_ 

We have discarded the second term in the integrand, which is exponentially subdominant 
except at X- = 0. For z > 2, the remaining integral is now dominated by values of X- 
near the one that maximises the function g{x_) = (1 — 5;_)[1 + x_{z — 1)]. This value 
is X* = {z — 2)/[2{z — 1)]; the fact that x* > for z > 2 justifies a posteriori why we 
were able to discard the second term from (32). Inserting x_ = x* into the integrand 
now gives 

lnG{z,N-^oo) = Nln(-^^^^Y z > 2, (35) 

up to subextensive contributions. The value x* has a natural interpretation: in the z- 
weighted ensemble it is the minimal value of the random numbers Xi in the box without 
the minima. The other numbers in this box are then distributed uniformly over [x*, 1]; 
setting z > 2 is (for A^ -^ cxd) sufficient to exclude any smaller values. In the box with 
the minima, on the other hand, random numbers from the whole interval [0, 1] occur, 
but the probability density is higher by a factor z on [0, x*] than on [x*, 1] (and uniform 
within these two intervals). 

We can use the expressions for G in the two regimes to compute 
the expectation value, (M), of the number of minima using (M) = 
Em MP{M, N)z^ / Zm P{M, N)z^ = dlnG/dlnz to obtain, in the limit A^ ^ oo, 

{M) = l ""! ,, (36) 



3.2. Mean field theory: Minima and Maxima 

Within the same mean field model, we can also compute the probability weights for 
configurations which contain maxima as well as minima. Call the number of minima 
Ml and the number of maxima M2, and the probability of such a configuration 
P{Mi, M2, N). The most direct way of obtaining this probability is from the counting 
argument outlined above. The result is 

f2N-Mi-M2-2\ (2N-Mi-M2-2\ 

I N-2 J + l iV-Mi-1 1+'^^^-^^^^-^ 
P(Mi, M2, A^) = 2 ^ ^ ^ ^ (37) 
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Ml ~ 1 \ 

M2 ~ iv\ 


Mi/iV = 1 
Ma/iV = 1 


Ml ~ 1 
M2 ~ 1 


Ml ~iV 
M2 ~ 1 

►- 



Z\ 



Figure 2. Sketch of the mean field phase diagram for the number of minima, Mi, and 
maxima, M^^ in the random energy landscape model. 



and can be explained as follows. Take the case where the Mi minima are in the left 
box; the prefactor 2 then accounts for the opposite case where they are in the right box. 
Now the M2 maxima can either be in the left or the right box. Suppose they are in the 
left box. Then in our construction of first drawing yi, . . . , y2N and then labelling them, 
yi,..., yMi and y2N-M2+i, ■ ■ ■ , y2N need to have label 'left' while i/Mi+i and y2N-M2 have 
label 'right'. The number of arranging the 2N — Mi — M2 — 2 labels that remain in the 
bag, of which N — 2 are 'right', is given by the first binomial coefficient in the square 
brackets in (37). The second term is constructed in the same way but with the labels of 
y2N-M2j • • • 5 y2N reversed: now A^ — Mi — 1 'left' and A^ — M2 — 1 'right' labels remain 
in the bag. This counting argument works while Mi < A^ — 1 and M2 < A^ — 1 (since 
we fix Ml + 1 labels at the bottom and M2 + 1 labels at the top). There is only one 
configuration that is not captured, namely. Mi = M2 = A^, where all labels are fixed: 
the third term of (37) accounts for this. 

For finite Mi and M2, where the third term of (37) is irrelevant, one easily sees 
that P{Mi,M2,N — > 00) = 2~*^i~^2- ^]-^g populations of minima and maxima are 
uncorrelated. (Configurations with the minima and maxima in the same and in different 
boxes also have the same weight, each contributing half the result.) Defining a generating 
function 



G'(^i,^2,iV)= J2 P{M,,M2,N)zf 



Ml M2 



Zn 



{31 



Mi,M2 



this implies G{zi, Z2, N -^ 00) = [-2i/(2 — zi)][z2/{2 — Z2)] for zi < 2 and Z2 < 2. In the 
ensemble weighted by zi and Z2 the average numbers of minima and maxima are then 
(Ml) = 2/(2 - zi) and (M2) = 2/(2 - Z2), respectively. 

For larger Zi or Z2 one can proceed using steepest descents as explained in Appendix 
A. We find that there is a first-order transition line at 1/zi + I/Z2 = 1. Beyond this {i. e. 
for 1/zi + l/z2 < 1), the numbers in the two boxes separate essentially completely, with 
one containing only minima and the other maxima: {Mi) /N = {M2) /N = 1. Between 



Phase Transition in a Random Minima Model 13 

this transition and the other boundaries aX zi = 2 and Z2 = 2 he two regions where the 
number of minima is extensive but the number of maxima is not, and vice versa. E.g. 
when Zi > 2 and Z2 < Z\j(^z\ — 1) one finds 

(Ml) = N (1^) , (39) 

IM\- ^^^' + ^^'^ ~ 2^12:2(2:1 + Z2) + zjzl 

^ '^ (2-Z2)(Zl-Z2)(^l + ^2-^1^2) ' ^ ^ 

with an analogous resuh when the roles of z^ and z^ are swapped. The last factor in 
the denominator for (M2) diverges at Xj z^ + 1/^2 = 1, signalling the transition to the 
regime where both Mi and M2 are extensive. Figure 2 shows a sketch of the overall 
phase diagram. The first order transition at Xjz^ + 1/2:2 = 1 has unusual features 
(see Appendix A): at the transition, an entire one-parameter family of phases becomes 
degenerate to leading order, i.e. has the same value of iV~^lnG'. This should produce 
unusual finite-size scaling effects which we have not yet explored. 

3.3. Bethe lattice 

Returning to the random minima problem, let us summarize the results so far. Within a 
two-box mean field theory, we found that the typical densities of minima in both boxes, 
p = {M)/N and p = {M)/N, vanish for fugacities 2: < ^c = 2. For higher fugacities, 
the system orders, with a nonzero density of minima p = {z — 2)/{z — 1) in one box 
but a vanishing one in the other, p = 0. A peculiar aspect of this behaviour is that the 
nonzero minima density increases linearly with 2: — 2:c = 2: — 2 around the transition, 
suggesting an order parameter exponent (3 = 1. For a mean field system this would 
be very unusual indeed as one would naively expect /? = 1/2 in mean field theory. We 
therefore next consider the minima problem on a Bethe lattice of finite connectivity. 
While the large connectivity limit should then retrieve the mean field results, at finite 
connectivity we would hope that a standard mean field phase transition with (3 = 1/2 
will reappear. This is indeed what we find. 

We begin as in the hard particle scenario by considering a Cayley tree. The basic 
quantity of interest is P^^\M^''\ m^''\ x^''^) , the probability - under random sampling of 
the Xi - that the root node of a Cayley tree of depth / has function value x^^\ that it is 
(or is not) a minimum as indicated by m^''^ = 1 (m*^'^ = 0), and that there are a total 
number M^''^ of minima in the tree. The basic recursion for this is 

X '^M('+l),...^m('+i),... (41) 

which expresses the fact that x^^^^' at the new root node is chosen independently of 
what happens in the p different branches i = 1, . . . , p attached to it. Once a;*^'"*"^^ and 
the properties of these branches are known, the values M*^'"*"^) and m^'-'^^^ for the new 
(/ + l)-level tree are fully determined as indicated schematically by the delta-functions. 
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Explicitly, 



m 



C+i) 



if xC+i) < 
otherwise 



X. 



(0 



M(i+^) = M« + m^+i) - ^mf 0(x 



Vi 



,/x 



„(0 



X 



a+i) 



(42) 



(43) 



The last sum runs over the /i nodes below the new root node as before. It expresses 
the fact that even if these nodes were minima within their own subtrees, once they are 
connected to the new root node they cease to be minima if they have function values 



^,(0 ^ T,(/+i) 



Introducing the generating functions for P*^'\ 



G«C 



X] 



M=0 



M 



p^^\m, 



m, X 



(44) 



where the fugacity z again acts on the number of minima, the recursion becomes (the 
restriction < a; < 1 is understood for all x- variables and so in particular Qix^^^^^') = 1): 



G^^i 



X] 



.(0 



dy G\;>{y)+G'l\y) + dy 



(0/ 



Gi;\y) + z-'G^I\y) 



dy G'^\y) + z-'G'l^y) 



(0/ 



G'l 



C+i), 



X] 



z dy G\l\y) + z~'G':\y) 



(0/ 



(45) 
(46) 



The second of these is easiest to explain: if the root node is a minimum with function 
value X, all fu, nodes in the level below must have function values x\ = y > x. The 
factor z~^ in front of G[ (y) corresponds to the negative term in (43), i.e. the fact that 
none of these nodes can then be minima. The prefactor z accounts for the new minimum 
at the root. The recursion (45) works similarly: multiplying out the /i-th power in the 
first line and subtracting the term in the second line gives all the possible configurations 
where at least one of the nodes below the new root has a lower function value than the 
latter. The factor of z~'^ is again for nodes which have higher function values than the 
new root and so cease to be minima if that is what they previously were. 

As in the hard particle case, a simpler recursion is obtained by taking ratios of 
appropriate generating functions. Here, it turns out to be convenient to consider the 
ratio S^\x) = G^!^ (x) / Gq (1) . We also abbreviate 



H^^\x)= fdy\s!,'\y) + z-'S?iy) 
/«(x) = rdy[S^\y) + Sf\y) 

Jo L 

Then our recursions read simply: 

{\^^Yst^\x) = i'^^\xY - h'^^\xY 



M) 



+ dy S^"{y) + z-'S['>{y) 



-1q(1), 



(x^'YSi 



(m), 



X) 



zH^'\ 



X] 



(47) 
(48) 

(49) 
(50) 
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where (X^''')^ = /'■'•'(l)^ is the normahzing coefficient that enforces 5*0 (1) = 1 for all 
/ as it must be. The corresponding differential versions will be more useful for later: 
expressing the S^^^\x) as derivatives of H^''~^^\x) and I'^^^^\x) gives 

(A('))^9,/J('+^)(^) = -I^^\xy (51) 

(A('))^9,/('+^)(a:) = {z - l)H^^\x)'' (52) 

with boundary conditions 

H^'\l) = 0, /j(')(0)=/(')(0) . (53) 

Once we have the functions H^''\x) and /*■'-' (x), the density {i.e. the probability of having 
a minimum) at the central node of a Bethe lattice follows directly as 

,a+i) = zjdxH(^\xr^^ 

^ jdx [j(o {xY+^ - m) {xY+^] + zjdx ij« {xY+^ ■ ^ ^ 

We expect as in the hard particle case that iteration of the above recursion over 
/ either gives /-independent values (disordered phase) or alternating layers (ordered 
phase): H^'^^\x) -^ H{x), if(^^+^)(a;) -* Six), and similarly for I^'-\x), A^'^ and p^^\ It 
is convenient to study directly the ordered case since it includes the other. The fixed 
point equations that one needs to solve are then 

X^d^H{x) = - I{xY (55) 

\^dj{x) ={z- l)H{xY (56) 

Vd^H{x) = - I{xY (57) 

X^dJix) ={z- l)H{xY (58) 

with the boundary conditions (53) holding for the functions in both the even and odd 
layers, and with 

A = /(I), A = 7(1) . (59) 

A very useful property that follows by combining the fixed point conditions is 

{z - l)X^H{xY+^ + X^I{xY+^ = A^A^+^ (60) 

independently of x: the x-derivative of the l.h.s. vanishes, and the value on the r.h.s. 
can be obtained by setting x = 1. This identity allows one to decouple the fixed point 
conditions for H and /, giving for the former 

d,H = -[l-{z- 1)X'^X''H^+Y^''^'^ ■ (61) 

Integrating by separation of variables and using the boundary condition H{1) = yields 
the following implicit expression for H{x) 

with B(p,q;a) = ^q dtt'P^^{l — tY~^ the incomplete Beta function. For H{x) one has 
the analogous result with A and A swapped. 
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It now remains to find A and A. To this end one can exploit the remaining 
conditions H{0) = /(O), H{0) = 7(0), from (53). Combining with (60) at x = 
and the corresponding relation with even and odd layers swapped, we find 



)M+i = i(oY+^ = y+ii^_zmN2iLil 



H{0)^+' = 1(0)'^+^ = A'^^^^ f-^ (63) 

and similarly for the odd layers. Inserting back into (62) for a; = gives 

, /2-l\ 1/(^+1) /I 1 , ,2-l-A/A\ , , 

The same relation again also holds with even and odd layers swapped. Together, these 
two conditions determine A and A. In the disordered phase, where A = A, A can be 
trivially found from (64) in closed form (and is equal to A(0) as defined below). 

For the densities, equation (54) suggests that one might need the explicit forms 
of H{x), I{x), H{x) and I{x). However, after some algebra one gets, by transforming 
integrals over x to integrals over H using (61) and similarly for integrals involving /, 
expressions that depend only on A and A, and indeed only on their log-ratio r = ln(A/A): 

2 1 . /^ T\z-l-e- 



-^B i^,^:(z-l 



5(^,^;(^-l)^^:^)+e-2'-/(M+i)E(^,^;(^-l)^fi^) ^ ^ 

V/i+l' /i+l' V / z(z-2) J \M+1' M+1 ^ ' z{z~2) J 

with an analogous expression for p. In the disordered phase, where r = 0, the single 
density is then given by the relatively simple equation of state 

^^^Uil^m^_ (66) 

'^ z-l B^ ^■^) 
^ Vm+1' M+1' 2 / 

To understand the solutions for A, A in the ordered phase, it is useful to have a 
single condition for r. Equation (64) gives A = A(r) with 

A(r) = (^ + 1) [{z - l)et/^^+^^ B-' (^, -±-; {z - 1)^^^) • (67) 

\/i+ 1 ;U+ 1 z{z - 2) ) 

The swapped relation gives A = A(— r) or A = e''A(— r). Since the two expressions for A 
have to agree, the desired condition on r is 

e-''^A{r) - e"/2A(-r) = . (68) 

The disordered phase has r = 0, which is the trivial solution. The bifurcation to 
the ordered phase takes place when the first r-derivative at r = vanishes, i.e. when 
(29rA — A)|^^Q = 0. This gives the condition 

\p+l p + l z^ j fi-l Zc-2 ^ ^ 

for the critical value z^ of the fugacity. It is easy to see that r initially departs from 
as (2; — ZcY^'^ as z is increased to above z^, this follows because (68) is odd in r so 
when the first derivative vanishes the leading term is third order in r. The densities (65) 
then have the same leading-order square root singularity. At generic finite connectivity 
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/i + 1 (> 2) we therefore retrieve, as hoped, an ordering phase transition with a standard 
mean field order parameter exponent [3 = 1/2. 

For generic z and fi one needs to solve numerically for r from (68) and then calculate 
the densities p and p from (65) and its analogue with even and odd layers swapped, i.e. 
r ^ —r. Further analytical progress can again be made for large p, however. We will 
need in particular the scaling of z^ for large p. One uses that 

for p ^ oo and a = 0{1) fixed (see (B.7)). The r.h.s. of (69) must then also diverge 
for p -^ oo, hence Zc -^ 2. To leading order the l.h.s. is /x + 0{1) while the r.h.s. is 
4/ (-2c — 2); this forces 

z, = 2 + A/p + Oil/p^) (71) 

for large p. 

3.3.1. Large p, above the transition We need to find first how r scales for fixed z > z^ 
and p -^ oo. Let us take r > for definiteness since solutions come in pairs (r, — r), 
and write (68) as 



B^.^-Az~l 



. z-l-e' 



gr(l-^)/(/.+l) ^ Vm+1'm+1'^ ' z{z-2) / ,^2) 



,^+1' ^+1' ^ ' z{z-2) 

For large p, the l.h.s. becomes e~^\ for this to be < 1, the third arguments of the Beta 
functions on the right cannot stay bounded away from or 1 since otherwise their ratio 
would converge to unity from (70). The third argument of the numerator Beta function 
thus has to approach zero, i.e. e"^ = z — 1 — 5r with 5r -^ 0. To leading order we have 
then 



z-l B -^ ^-1 ^-^ 

-^ Vm+1' M+1' (z-1)z{z-2) 



^ Vm+1' /^+1' 2(^-2) / /-jz^N 



i?(^,^ -5(^,^ 



5r 



Now for a remaining finite or going to zero, (/i + 1)^^-8 f^j-) ^i! '^j ~^ a^l'^^^^^ 
for p ^ oo (see (B.4)). Bearing in mind that the complete Beta function obeys 
{p + l)^^B (— h", -^) -^ 2, we get to leading order 



I (5^1/(^+1) 



^ _ 1 2 - 5rV(M+i) C^S) 

or, calling a the limiting value of Sr^'^^^^', a = 2/z; 5r thus decays exponentially with 
p as 6r ^ {2/ zY. 

The density in the even layers can now be worked out from (65). To leading order, 
using that 5 (gf , ^; 1 - a) = 5 (gf , ^) - B (^, gf ; a)^(p + 1)(1 - aV(.+i)) 
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(see after (B.4)) 



^q(b±1 ^. 1 
2-1-^ Vm+1' m+1' 



5r 



{z-l)z{z-2) 



-^ Vm+1' m+1' (2-1)2(2-2) y ^ VM+i' M+i' 2(z-2) 



2; 1 — a z — 2 



(76) 
(77) 



z — 11 — a + a 2;— 1 

This agrees with the simple two-box mean field theory as we had hoped. 

The density in the odd layers, on the other hand, goes to zero for large /x: 

_z_Tj(lJ±2 1 . {z-l)Sr \ 
^ ^ z-l Vm+1' M+1' ziz-2) J .„ 



R fA£±2 ^_. (2-l)^r \ r^ /_|_ £+2. -, 5r 

^ Vm+1' m+1' 2(2-2) ) '^ VM+1' M+1' (2-1)2(2-2) 

5r 



(79) 



(2;-2)(^ + l) • 
This is exponentially small in /x because 5r is. 

We discuss briefiy the behaviour for large /x of 5*0(3;) and Si{x). Up to an overall 
normalization factor these give the probabilities, in the random landscape ensemble 
weighted by the fugacity z, that the root node of a Cayley tree has random function 
value X and is (for 5*1(0;)) or is not (for 5*0 (x)) a minimum. For large /i, the fact that 
the root node has /i rather than ^ + 1 neighbours becomes unimportant and these 
probabilities also apply to an arbitrary node in the bulk of the Bethe lattice. 

We find by constructing the explicit solutions for H{x) etc. that there is a threshold 
value of X, x* = {z — 2)/[2{z — 1)] so that for yU ^ 00 one has 5*0(3;) = 5*0(0;) = Q{x — x*), 
5*1 (x) = zQ{x* — x), 5*1 (a;) = 0. (Correspondingly, the functions H and / are piecewise 
linear below and above x*.) So in the even layers, i.e. those with a nonzero density of 
minima, a site is a minimum if its value Xi is below x*, and not a minimum otherwise. 
In the odd layers, no sites are minima, and the values Xi at all sites are above x*. This is 
exactly the same phenomenology as in the two-box mean field theory, confirming again 
that the latter captures most of the physics of the large connectivity limit on the Bethe 
lattice. The exception is the region around the ordering transition where the square-root 
singularities and hence the order parameter exponent (3 = 1/2 are visible: this becomes 
vanishingly small as we will now see. 

3.3.2. Large fi, around the transition From the large-/i expansion (71) we expect that 

the appropriate scaling for the fugacity in the region around the phase transition is 

z = 2 + z/fi. As we will see, this corresponds to r being of order l//i, r = r/fi. With 

these scalings, the third arguments of the Beta functions in (72) become |(1 =F r/z). 

Using from (B.7) that for /i -* 00, 5 (^, ^; |(1 + a)) - (/i + 1) ^ ln[(l + a)/(l -a)] 

and keeping only terms of (9(1) and 0{l/fi) gives 

f /i + 1 + lnf(l -f/5)/(l +f/5)l 2 , ,, , , ,, 

1 - - = ^ . ' {';, ' = 1 + - n 1 - f /5 / 1 + f/5 . 80 

Equating the (9(l//i) terms shows 

f = 2 ln[(l + f/z)/{l - f/z)] = 4 artanh(f/5) (81) 
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Figure 3. Equation of state of the random minima problem around the ordering 
transition, for the Bethe lattice in the limit of large connectivity /i + 1. Shown are the 
scaled densities of minima, /?, in the two boxes, against the scaled fugacity z. Inset: 
After transforming nonlinearly to m = 1 — 2e~^/^, the equation of state becomes that 
of a mean field ferromagnet at inverse temperature 5/4. 



or 



2) 



tanh(r/4) 

The critical point is reached for f ^ 0, giving ic = 4 in agreement with (71). 

It remains to work out the densities. In the expression (65) for the even layers, the 
arguments of the Beta functions simplify as before, and also zji^z — 1) ^ 2, so that 



9R f££±2 ^_. 1 



B 



i£±2 ^_. 1 / 1 



+ 5 



1 M+2 . 1 
M+1' M+1' 2 



1-? 



(83) 



The second Beta function in the numerator equals /i + 1 to leading order while the other 
ones are 0(1), B (^^, -r^-, a) — ^ /o"'^^ (1 ^^) ""^ = — ln(l — a), so that the scaled density 
p = Pli becomes for yU ^ oo 



P 



-2 In 



2 V z 



and similarly in the odd layers, after swapping r -^ — r, 
y^ = -2 In 



2 V 5 



^4) 



^5) 



The equations (82,84,85) give the equation of state in the phase transition region. The 
occurrence of the auxiliary parameter f is a little awkward but can be eliminated if we 
transform the (scaled) densities nonlinearly as 

- 2e-~p''^, u=l- 2e-^/2 (gg^ 



u 
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so that, from (84,85), u = —u always. By combining (82,84) one then sees that 

u=^- = tanh(f/4) (87) 

and so finally z = r/u = 4artanh(-u)/-u or 

u = tanh((i/4)M) . (88) 

Since u = —u, the same equation also holds for the density in the odd layers. 
Remarkably, therefore, once the densities are nonlinearly transformed according to (86), 
they depend on the fugacity exactly as the magnetizations in a mean field ferromagnet 
with unit interaction strength and inverse temperature z/4. We show the equation of 
state in Fig. 3, both in terms of the (scaled) densities p and, in the inset, the transformed 
variables u. 

3-4- Mean field theory with soft minima 

In this final subsection we ask whether the behaviour around the ordering transition 
that we found for a highly connected Bethe lattice can also be obtained directly within 
a mean field theory. It turns out that this is possible: drawing inspiration from our 
treatment of the hard particle model, we make the labelling of sites as minima 'soft'. 

The two-box setup is initially the same as for hard particles, with generating 
function 

G{z, N) = (z^+*^) = (^z^M+^^)) (89) 

Here the average is over our random landscape ensemble as before, while M and M label 
the total number of minima in the left and right box, respectively. For hard minima 
of course one and only one of these quantities is ever nonzero; for soft minima both M 
and M can be nonzero. 

To define 'soft' minima, we first introduce an auxiliary variable Tj G {0, 1} at each 
site. We can obtain the usual generating function for hard minima by forcing this to be 
if mj = 0; otherwise we allow it to be or 1. Assigning weight factors 1 and z — 1 to 
Tj = and 1, respectively, we can then write the factor from each site in the generating 
function as 

z-^ = E (^ - l)"'^.,(i-™,),o . (90) 

n=0,l 

(Indeed, for mj = only Tj = is allowed and we get z^ = 1 = {z — 1)°; in the opposite 
case we have z^ = {z — 1)^ + {z — lY.) Now to make the minima soft, we relax the 
constraint that Tj = if m^ = 0, i.e. we replace 6r^(^i-m^),o — ^ exp[— arj(l — rhi)]. Here 

N 

mi = N-^Y.^(^3-^i) (91) 

i=i 

is a soft version of mi. it measures what fraction of numbers in the other box are above 

Xi, so that rhi = 0, 1/A^, . . . , 1 — 1/A^ corresponds to rui = and rhi = 1 to rui = 1- 

Thus, when rrij = we have 1 — rhi > 1/-^ and for a — ^ oo at fixed A^ our soft minima 
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weight exp[—aTi{l — rhi)] reverts to 5T--(i_mi),o ^^ it should. As in the hard sphere case 
we in fact take N —^ oo first and then a -^ oo. 

We summarize our starting point: the generating function for soft minima is 

G{z, N) = Tr,,^ {z - l)Eti(-»+->)e^^ (92) 



A = — In ( exp —a T^ 



riil-mi) + fiil-ifii)]\) (93) 



where Tr^-.^ abbreviates the sum over all Tj and fj. The soft version of the minimum 
indicator variables rhi in the right box is defined in the obvious way by swapping the 
roles of X and x in (91), i.e. rhi = N~^ Y^jLi Q{xj — Xi). 

To calculate G, consider first the average from (93). By permutation symmetry 
within each box, this can only depend on the numbers T, T of nonzero r's in the two 
boxes. Writing out the definition of the rhi and rhi in terms of sign functions, this gives 




a 



2N 



T N T N 

im(l + sgn(a;i - Xj)) + ^ ^(1 - sgn(a;i - Xj)) 

1=1 j=l j=l i=l 



(94) 



= -a^^ + biT,N-T) + biN-T,T) (95) 

Here we have used that the sgn terms with 1 < j < T in the first sum exactly cancel 
those with 1 < i < T in the second one, so that the remaining average factorizes into 
two independent terms of the form 

e^^(^-^^) = ^exp (-^EE ^Sn{x.. - ^.)] ) • (96) 

The replacement Xi ^ 1 — Xi, Xj ^ 1 — Xj leaves the distribution of these variables 
unchanged, hence b is symmetric under a — > —a, as well as under interchange of Ti and 
T2. To evaluate b, we can assume without loss of generality that the Xi are ordered. 
Using also that the average over the Xj factorizes. 



giVfe(Ti,T2) 



+ (l-XTje^^"/2^]'^') (97) 

where the remaining average is over the Xi. Denote the quantity raised to the power T2 
by y. Setting also Vq = Xi, Vi = 0:2 — Xi, . . . , Wti-i = ^t^ — Xtj^-i, f Ti = 1 — Xti, the 
Vi are non-negative (because of the ordering of the Xi) and uniformly distributed apart 
from the constraint J2i=o Vi = 1, such that -P({f «}) = Ti\ 5{l — J2i ^i)- The characteristic 
function of y is then 

^Nu.y\ =T^i I !1^ I^Y[dv,expi N\{1 - V t;,) + iVcc; V t;,e(2-^i)"/2^ 1 (98) 



J Zm JO j=o \ j i 

T,\ f ^e^" n (^^ - iVcue^*-^^/^)"/^)"' . (99) 



Phase Transition in a Random Minima Model 22 



Reverse Fourier transforming now produces 

Ndiu r NdX 



^Nb{TuT2) ^ /yT2 



T,\ fdy I -^ f cxp 

J J 27ci J 27ri 



T2\ny-Nujy + NX 



Ti 



- (Ti + l)lnA^-^ln(A 



cue 



{i~Ti/2)a/N 



(100) 



j=0 

Defining the intensive quantitites ti = Ti/N, t2 = T2/N, we can do the integral using 
steepest descents ior N ^ 00: 

b{ti,t2) =max\ti\n{ti/e)+t2lny-uy + \- f rfw ln(A - cue"") I . (101) 

y,'^A [ J-ti/2 J 

Setting the derivatives w.r.t. y, u and A to zero gives the saddle point equations 

L0 = -, y= du- , 1= du- . (102) 

y J-h/2 X-uje""" J-ti/2 A-c^e"" ^ ^ 

Combining the last two we find 

,.ir,„(^^_i).A^ (103) 

uj J-ti/2 \A — cue"" / u 

and hence A = ti + t2. In the last saddle point equation we can perform the integral 
explicitly, yielding 

^^±,J^'-"f--] (104) 



and we can solve for u: 



3^*2/2 _ p-at2/2 



Together with \ = ti+t2, y = t2/uj we thus have all saddle point values explicitly. The 
derivatives of b{ti,t2) that we will need become 

^ = Inti - - In (A - cje"*i/2)(^ _ ^e'^^^/^) ^^06) 

Oil 2 L J 



lT^ + ^+l"l .\-„., (107) 



ti at2 , /l_e-"(*i+*2)^ 

1 h In 

ti + t2 2 \ l-e- 

and db/dt2 has the same form with ti and t2 interchanged. This symmetry 
property is clear from the definition (96) but not so obvious from the saddle point 
representation (101). An explicitly symmetric expression can be obtained from (107) 
by integrating from ti = (where 6 = 0); after a little algebra, this can be cast in the 
form 

ab{tut2) = F{a{ti + ts)) - F{ati) - F{at2) (108) 

with F{x) = Jq du ln[2sinh(-u/2)/-u]. One might hope that a derivation exists which 
directly reveals this simple structure, but so far we have been unable to find one. 
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Now we can finally write down the saddle point equations for the full generating 
function. In terms of the r-densities t = T/N, i = T/N, one has from (92,93,95) 

^lnG= max{n{t)+n{t) + {t + t)\n{z-l) + A} (109) 

= max{nit)+nii) + {t + i)[\n{z- I) - a/2] + b{t,l - i) + b{l - t,i)} (110) 

where 7i (t) and 7i (t) again account for the combinatorial (entropic) contributions. 
Then 

db 



9 1,^ , 1 — t , , _, , ,^ db 

__,„G=,„_ + ,„(,-l)-„/2+_ 



ti=t,t2 = l-t "^^l 



:iii) 



ti=l^t,t2=t 

l-t + i\ , /I - e-"(i+*-*^ 1 - e-"(i"*) ~ 



= ln(z-l)-at + ln ^ + In ——j- — (112) 

and this must vanish at the saddle point. The corresponding equation for t just has t 
and i swapped. 

In the disordered phase t = i, the saddle point equation can be solved explicitly to 
get t = a~^ ln[2;/(l+e~")]. One can then again ask about bifurcations to solutions where 
t ^ i. The critical value of z can be got as follows: think of the first saddle point equation 
as defining implicitly t as a function of t; the second saddle point equation defines the 
inverse function, which graphically is flipped about the diagonal. The disordered flxed 
point on the diagonal becomes unstable when the slope di/dt = —1. The resulting 
condition on z^ looks complicated, but neglecting terms that are exponentially small in 
a one gets Zc = 2(1 — q;^^)^/(1 — Aa^-^) = 2 + 4a~^ + 0{a^'^). This looks encouraging: 
with the identiflcation fi = a, it is identical to the result (71), suggesting that the soft 
minima mean fleld theory captures the large connectivity limit on the Bethe lattice. 

We recall for the evaluation in the following subsections that t = N~^J2iTi^ ^^^ 
similarly i, are the densities of the r- variables. For large enough a, we can have r, = 1 
only when there is genuinely a minimum at site i {irti = 1); but even if m, = 1 then 
Tj = with probability 1/(1 + z — 1) = 1/z. So for a ^ oo the r-densities are related 
to the true densities of minima by t = [{z — l)/z]p, i = [{z — l)/z]p. 

3.4-1 ■ Large a, above the transition Here we expect that one of the two densities (say 
t) stays flnite and < 1 while the other (t) goes to zero. The saddle point equations are 
then, up to exponentially small terms: 

0= \n{z-l)-at + \n( ^'^^J^ (113) 

= ln(z - 1) - at - In ( 1^A±| J _ in(l _ e""*) . (114) 

In the second equation, the only way to balance the —at term is to have i vanish faster 
than 1/a so that the argument of the last log tends to zero; the log itself can then be 
approximated as — ln(Q;t). This gives to leading order i ~ Q;~^e""*. Inserting into the 
flrst equation then leads to (l + t)/(l — t) = z — 1 or t = (z — 2)/z. The density of minima 
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in this box is therefore p = [z/{z — l)]t = {z — '2)/{z — 1), consistent with our direct 
calculation in the hard minima limit. The critical fugacity is Zc = 2, also as expected. 
Overall the a ^ oo limit correctly reproduces the hard minima scenario as desired. 

3.4-2. Large a, around the transition By analogy with the Bethe lattice calculation, we 
scale the r-densities and the fugacity as t = i/a, i = i/a and z = 2 + z/a, respectively. 
The saddle point equations are then, again up to exponentially small terms, 

'l-{i-F)/a, , ,^ , 



0= In(z-l) -t + ln ( ^ — J^— I -ln(l-e"*) (115) 

1 + (t — i)/a ^ 



0= ln(z-l)-t-ln ( V J' I -ln(l-e"*) (116) 

1 + (t — t)/a 

It is again useful to make a nonlinear transformation from the r-densities to 

M=l-2e-*, M = l-2e-*. (117) 

The first saddle point equation then implicitly defines a function U via u = U{u). The 
second one gives u = U{u). For these to be consistent with each other, we require 

a[U{u)-U-^{u)] = {) . (118) 

Here U^^ is the inverse function of t/; the factor a will be useful shortly. To find the 
function U for large a, we expand the first saddle point equation, keeping terms of 0{1) 
and 0{l/a): 

z ^ 2{i-~i) , ^^ _f^ ^ fl-u\ ^ z + 2\n[{l-u)/{l-u)] 



= t--^ ^-In l-e^* =ln + ^ '-^ ^. (119) 

a a ^ ^Vl + M/ a 

To leading order this gives u = U{u) = —u: the function U is identical to its inverse. 

This is why we need to go to 0{l/a) to get a nontrivial condition for u, as emphasized 

by the factor a in (118). Now insert the leading order relation u = —u into the 0{1/q) 

term above to get 

,-. = f/M = i-(i + ,.)(i-i±lMi_^Mi±i^) (120) 

a 
The inverse function is obtained by solving (119) for u: 

,„ . [,-!(,) =_!+(!_ ,,) (i + S + 2M(l + u)/(l-.)lj i^^^j 

= — -u + (1 — u) . (123) 

a 

So the equation (118) determining u becomes for a ^ oo 

= (1 + m)(5 + 2 ln[(l - u)/{l + u)])-{l-u){~z + 2 ln[(l + u)/{l - u)]) (124) 

= 2zu + A ln[(l - m)/(1 + u)] =2zu-% artanh(M) (125) 
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or 

u = tanh((5/4)M) . (126) 

This is exactly as on the Bethe lattice around the transition, so the entire scaling 
behaviour in this region matches between the two cases, namely, the Bethe lattice in 
the limit of large connectivity /i and the soft minima problem in the nearly hard limit of 
large a. Notice that, while the definitions of the relevant nonlinear transformations (86) 
and (117) of the density variables look different, they are in fact identical because 
p= [z/{z-l)]i=2i+0{l/a). 

4. Summary and outlook 

In summary, we have analysed the number and distribution of minima in random 
landscapes defined on non-Euclidean lattices. Using an ensemble where random 
landscapes are reweighted by a fugacity factor z'^ depending on the number of minima 
M, the simplest viable (two-box) mean field theory showed an ordering phase transition 
aX Zc = 2. For z > Zc, one box contains an extensive number of minima with density 
p = {z — 2)/{z — 1). The onset of order seemed to be governed by an unusual order 
parameter exponent /3 = 1, which motivated our study on the Bethe lattice. 

Using recursion techniques, we found a full solution of the problem on the Bethe 
lattice which showed that for any finite connectivity fi + 1 (> 2) there is indeed an 
ordering transition with a conventional mean field order parameter exponent f3 = 1/2. 
As p becomes large, the region around the transition where this behaviour is visible 
shrinks as 1/p. It disappears as /i -^ oo at fixed fugacity z, and this is what causes the 
unusual effective exponent in the two-box mean field theory. We analysed separately 
the scaling for large fi for fixed z above the transition and for z within 1/fi of z^. In the 
latter case, a nonlinear transformation turns out to map the equation of state neatly onto 
that of a mean field ferromagnet. Finally, we showed that the region around the phase 
transition can also be analysed directly within a mean field approach, by making the 
assignment of minima 'soft' and then taking the nearly hard limit {a -^ oo). This was 
motivated by our analogous treatment of the hard sphere lattice gas, where a softening 
of the nearest neighbour exclusion revealed the ordering phase transition that remains 
entirely hidden within the two-box mean field theory. 

In the mean field approach we also considered the joint distribution of minima and 
maxima of random landscapes. Here two fugacities enter, zi and Z2, and in addition 
to the phase transitions ai zi = 2 and Z2 = 2 where the number of minima and 
maxima respectively first becomes extensive, there is a first-order transition on the 
line 1/zi + l/z2 = 1. Beyond this line, essentially all points in the landscape are either 
minima and maxima; in our mean field setup, these sites are separated into the two 
boxes. 

In future work, it should be possible to extend the analysis of joint distributions 
of minima and maxima to the Bethe lattice. This would presumably require three 
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generating functions, for sites that are minima, maxima or neither. Generahzing the 
soft minima/maxima approach looks less easy because for 'soft' labels one would also 
have to consider sites that are labelled as both minima and maxima. It would also be 
interesting to generalize further, and consider not just minima but also nodes with fixed 
number /c = 1, 2, ... of lower-lying neighbours. 

Finally, one would like to extend our calculation also to large random graphs with 
the same local structure as a Bethe lattice, i.e. regular graphs where all nodes have the 
same number (/i + 1) of neighbours. Given that short loops are rare on such graphs, one 
might intuitively expect to see the same phenomenology. However, the strict sublattice 
ordering on the Bethe lattice cannot be maintained in the inevitable presence of at least 
some loops with an odd number of links, and so in actual fact it is likely that one would 
instead obtain glassy phases as in related hard particle models [34, 35]. Generalizing 
our approach to this scenario appears to be a challenging problem indeed. 
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Appendix A. Minima and maxima for zi > 2 or Z2 > 2 

We outline two methods for understanding the two-box problem in the case where we 
track both minima and maxima. The first one starts from the large-deviation form 
of (37). It is easy to see that the first and third terms can never be larger than the 
second. Taking A^ large at fixed densities pi = Mi/N and p2 = M2/N of the minima 
and maxima then gives up an irrelevant constant 

N-HnP{Np^,Np2,N) = {2 - p^- p2)n( ^ ^ ~ ^' ] . (A.l) 

\z - Pi - P2 J 

If one rewrites the definition (38) of the generating function as an integral over pi and 

P2, the latter will therefore be dominated by those values maximizing the function 

7(Pi, P2) = (2 - pi - p2)'H ^ — + pi\nzi+ p2 In Z2 ■ (A. 2) 

V2-P1-P2; 

Now take for definiteness zi > Z2, so that any maxima will obey pi > p2- In this regime 

we can set 1 — pi = k{2 — pi — P2) with < k < 1/2 and have at fixed k a /inear variation 

with 1 — P2: 

7(Pi, P2) = ln(2;i2;2) + (1 - P2)s{k) (A. 3) 

sU) = InUzi) -InUl - k)z2] . (A.4) 

I — K 

The slope function s{k,) now tells us where the maxima of 7(pi,p2) are. First we 

maximize over k; if the maximum value of s{k,) is positive, we get a maximum of 

li.Pii P2) at P2 = and hence pi = (1 — 2k) /{I — k), otherwise a maximum at p2 = 1 

and pi = 1. 
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The derivative of s(/t) is s'{k,) = —ln{K,zi)/{l — k,)"^. For zi < 2 this is always 
positive and the maximum is at /t = 1/2, where s(l/2) = — ln(ziZ2/4) > (given that 
Z2 < zi < 2). So 7(pi,p2) is maximal at pi = P2 = 0, consistent with the analysis in 
the main text that showed that in this regime minima and maxima are both intensive 
in number. 

For Zi > 2, the maximum of s(/t) is at /t = 1/zi, where s{l/zi) = — ln[(l — 
l/zi)/{l/z2)]. If l/z2 > 1 — l/zi, this value is positive and 7(pi,p2) has a maximum 
at p2 = 0, pi = (1 — 2k,)/{1 — k) = {zi — 2)/{zi — 1). This is the mixed regime, with 
Ml extensive and M2 intensive, see (40). For 1/ Z2 < 1 — l/-2i, finally, the maximum 
value of s{k) is negative and 7(pi,p2) has its maximum at pi = P2 = 1- This is 
the fully separated regime, where one box contains essentially only minima and the 
other only maxima. Note that, as stated in the main text, at the first-order transition 
1/zi + 1/^2 = 1, the entire line in the (pi,p2) plane corresponding to k = Ijzi is 
degenerate, i.e. has the same value of 7(pi,P2)- A further peculiarity is that there is 
no metastability: neither of the phases persists as a local maximum of 7(pi,p2) on the 
corresponding 'wrong' side of the transition line. 

It remains to find the average number of maxima in the mixed regime {zi > 2, 
1/zi + l/z2 > 1). We already know that pi is nonzero then; on general grounds its 
fluctuations (~ l/\/iV) niust become negligible for large A^. We can then take the limit 
A^ ^ 00 in (37) at flnite M2 and Mi = Npi to get for the distribution of M2 at given pi 

Multiplying by ^2 ^, normalizing and taking the average of M2 then gives the result 
stated in (40). Note that (A. 5) has a simple interpretation in the labelling picture: 
For the Mi = A^pi minima we have used up as many 'left' labels (and one 'right' 
label). The probability that the largest number y2N will be labelled 'left' is then 
(A^ - Npi)/{2N - Npi - 1) ^ (1 - pi)/(2 - pi) for large A^. The flrst term in (A.5) 
thus gives the probability that y2N,y2N~i, ■ ■ ■ ,y2N-M2+i are all labelled 'left' and the 
next number down, y2N-M2) is labelled 'right'; the second term gives the analogous 
contribution from the reverse labelling. The respective probabilities (1 — pi)/(2 — pi) 
and 1/(2 — pi) for a 'left' and 'right' label remain the same throughout as we are only 
labelling finitely many (M2 + 1) numbers and so for large A^ the fraction of labels in the 
bag of either kind only changes negligibly. 

The second approach parallels more closely the one taken for the minima problem. 
We first calculate the joint probability distribution, P{x-,x+) of the smallest and largest 
number in the left box. The probability that all the xi are greater that some value x_ and 
smaller than some other value x+isV = (a;+— a;„)^. This is also the probability that the 
minimum of these numbers is larger than x_, and the maximum smaller than a;+, so the 
joint probability density of the minimum and maximum is obtained by differentiation 
as 

P(x_,x+)= -^^^ = N{N-l){x+-x.Y-\ (A.6) 
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There is an analogous expression for P{x_, x^). The remaining N — 2 numbers in each 

box are then distributed uniformly between x_ and x^, and a;„ and a;+, respectively. 

One can now represent the probability of getting Mi minima and M2 maxima in 

terms of averages over this distribution. As before we assume that the minima are in the 

left box, i.e. a;„ < a;„, and multiply the probability by a factor 2 to cover the opposite 

case: 

1 / (N-2V 

-P{Mi,M2,N) = {Q{x_-x_)Q{x+-x^^ ^ ' 



2 ' ' ' ' \ ' ' '^ ^'(M2-l)!(iV-Mi-M2)!(Mi-l)! 




M2-I /- - \ N-M1-M2 /_ \ Mi-1\ 

_[_ tXj — \ / eXy tAj - — \ 



a;_)0(a;+ — a;_)0(a;+ — a;+) 



M2-I 



x+ 


— x^ 


x+ 


— x^ 


x+ 


-x+ 




^Mi,n5m2,n {Q{x~ — xJ)Q{xj^ — x^)Q{x^ — x^)) . (A. 7) 

The three terms on the r.h.s. are arranged in the same order as in (37), and represent 
different orderings of x_, a;_, x+ and a;+. In the first term, the Theta functions and 
the constraint X- < x+ enforce the ordering a;_ < a;_ < a;+ < a;+, so that the left box 
contains both minima and maxima and the right box neither. The remaining factors in 
this term give the probability that out of the N — 2 numbers in the left box (other than 
x^ and x^) exactly Mi — 1 are below a;_ and hence minima, and M2 — 1 are above x^ and 
therefore maxima. The second term corresponds to the ordering x^ < x^ < x^ < x^, 
where the maxima are in the right box but the ranges of numbers in the two boxes still 
overlap. In this case we need to find Mi — 1 numbers in the left box (in addition to x_) 
that are below a:_, and M2 — 1 numbers in the right box (in addition to x^) that are 
above Xj^. Finally, the last term is for the ordering X- < x+ < X- < x^. All numbers 
in the left box are then smaller than in the right one, and we have Mi = A^ minima on 
the left and M2 = A^ maxima on the right. 

In the representation (A. 7) one can easily perform the sums defining the generating 
function (38) to get 

-G{zi,Z2,N) = / 0(x_ - a;_)6(a;+ - x+) 



X Z1Z2 



2:2 (x+ — x^) + x^ — x^ + Zi{x- — x^)\ ^ 
x+ — x^ I , 



+ ( 0(a;_ — x)Q{xj^ — x_)Q{x^ — x^) 

I X^ ~ X_ + Zi{x_ — X_)\ f Z2{x+ — Xj^) + Xj^ — X^\ ^ 

\ X-|- — X— I \ X-f- — X— I j 
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+ z^z^ {e{x_-x_)Q{x+-x+)e{x_-x+)) . (A.8) 

To carry out the averages one inserts (A. 6) and its analogue for the right box and 
integrates over a;_, x_, x^ and x+. In each term two of the integrals can be done 
directly and one is left with 

-G{zi,Z2,N)= dx- dx+N{N-l){x+-x^)'^~'^\\z2{l-x+)+x+-x^ + zixJ^ 
2 Jo Js_ '- 

— [2:2(1 — x^) + x^ — a;_]^ — [x^ — X- + Zix^]^ + [x+ — X-]^ \ 

+ f dx- f_ dx+ N"^ {[x+ - x_ + 2ix_]^"^ - [x+ - S_]^"^} 

X {[22(1 - x+) +X+- x^f-^ - [x+ - x^f~^} 

+ z^z^ C dx+ t dx. N^x1-\l - x.f-^ . (A.9) 

JO Jx+ 

The remaining integrals in the last line can of course also be done and give [(2A^)!/A^!^]~^ 
as expected from (37). 

From here on one can proceed as in the minima-only case. If both Zi and Z2 are 
below 2, one rescales x_ = u/N , x+ = 1 — v/N in the first integral and similarly for 
the other terms; this gives G{zi^ ^2, N -^ 00) = [zi/{2 — zi)] [-22/(2 — Z2)] as derived by a 
different route in the main text. For larger zi or Z2 one uses steepest descents again. The 
functions to be maximized always have negative Hessian determinants so the maxima 
are on the boundary. We illustrate only the mixed case zi > Z2 at 1/2:1 + 1/2:2 > 1. Here 
the relevant saddle point in the first integral is x_ = {zi — 2)/[2(2:i — 1)] and x^ = 1. 
Because this is at the upper extreme of the integration range of x^, however, one needs 
to rescale x+ = 1 — v/N to treat the near-cancellation of the first and third terms in 
the integrand explicitly. (The second and fourth terms make exponentially subleading 
contributions.) In the second integral one needs to set similarly a;+ = 1 —v/N to capture 
the near-cancellation in the second factor. One thus gets for large A^, after neglecting 
the exponentially subdominant third term of (A.9): 

1 rl /-oo 

-G{zi,Z2,N)= dx. dv{N-l){l-x.f-^e-''/^^-^-^ 
2 Jo Jo 

X [1 + (zi - l)x_]^|e('^2-l)f/[l+(^i-l)5;_] _ g-^/[l+(2i-l)S_]| 
rl roo _ 

+ dx. dvN[l + {zi-l)x.f'h-''/^^+^'-'~^^^-^ 
Jo Jo 

X (1 - x.f-' |g(.2-iW{i-5=-) _ g-./{i-5=-)| (A.IO) 

The common exponential factor {(1 — a;_)[l + (2:1 — l)a;_]}^ means that for large A^ we 
can replace x_ = {zi — 2)/[2(2:i — 1)] in all other, slowly varying, terms to obtain 
1^/., ., ^n 2N{z^-lfz2iz2-2) /-i^^ ^^^ ^ ,^^ , ,_^ ,^. n^ 



-Giz,,Z2,N)= . ' '/ '' ' '— / dx.{{l-x.)[l + {z^-l)x.]y' .(A.ll) 

2 2:1(2:1 -2:2) (2:12:2 - 2:1 - 2:2) Jo 

This result is asymptotically exact for A^ ^ 00. It depends on 2:2 only through 

subexponential factors, which is why the cancellations referred to above have to be 

treated so carefully. Using (M2) = dlnG/dlnz2 then retrieves after a little algebra the 

result for the number of maxima in the mixed phase stated in (40). 
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Appendix B. Beta function asymptotics 

Here we gather the asymptotic properties of the incomplete Beta function B{p, q; a) = 
jQdttP~^{l — ty^^ that we need in Sec. 3.3. Specifically, setting p = l/(/i + 1), we 
require the behaviour of B{p,p; a), B{p,p + 1; a) and B{p + l,p; a) in the limit p ^ 0. 
Directly from the definition one sees that these three functions are linked by the simple 
sum rule 

B{p,p;a) = B{p,p + 1-a) + B{p + l,p;a) . (B.l) 

We will always keep a < 1, with a either fixed as p ^ or itself going to zero. 
The last function in (B.l) is simplest as it remains non-singular: 

B{p + 1, p; a) = r dt t^i,! - tf-^ -^ /" rft (1 - t)-^ = - ln(l - a) . (B.2) 
JO Jo 

The other two functions, on the other hand, diverge as, p -^ 0. With the variable 

transformation s = t^ one gets 

pB{p,p;a) =p rdttP~\l-ty-^ = r ds (l - s^/py^ . (B.3) 

Now if a vanishes quickly enough (exponentially in 1/p) when p ^ for a^ to stay 
bounded below 1, we can exploit the fact that the integrand approaches unity for all 
s < 1 to get 

pB{p,p;a) ^a^ . (B.4) 

From (B.l) and (B.2), which shows that B{p+l,p; a) stays finite, the same limit applies 
to B{p,p + 1; a). 

The result (B.4) does in fact extend also to a that vanish more slowly or stay finite, 
so that a^ — !► 1. One can see this by subtracting off the leading term: 

B(p,p;a)--aP= f dtf^^^Hl -tf-^ -I] . (B.5) 

p Jo 

The integrand now remains non-singular for p ^ and approaches 1/(1 — t), so that 

B(p,p;a)--aP ^ -\n(l- a) . (B.6) 

p 

Multiplying by p gives back (B.4) as claimed. For fixed a it is more useful to rewrite 
the last relation, using a^ = 1 -|- p In a -|- O(p^), as 

B(p,p;a) >\n\a/(l-a)] . (B.7) 

p 
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